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Abstract. We propose and test a new method based on Richardson-Lucy deconvolution to reconstruct three-dimensional gas 
density and temperature distributions in galaxy clusters from combined X-ray and thermal Sunyaev-Zebdovich observations. 
Clusters are assumed to be axially symmetric and arbitrarily inclined with respect to the line-of-sight. No equilibrium assump¬ 
tion other than local thermal equilibrium is needed. We test the algorithm with synthetic observations of analytically modeled 
and numerically simulated galaxy clusters and discuss the quality of the density and temperature reconstructions in idealised 
situations and in presence of observational noise, deviations from axial symmetry and cluster substructure. We find that an¬ 
alytic and numerical gas density and temperature distributions can be accurately reconstructed in three dimensions, even if 
observational noise is present. We also discuss methods for determining the inclination angle from data and show that it can 
be constrained using X-ray temperature maps. For a realistic cluster and including observational noise the three-dimensional 
reconstructions reach a level of accuracy of about 15%. 


1. Introduction 


In hierarchical models of structure formation, galaxy clusters 
are not only the most massive gravitationally bound objects in 
the Universe, but also the most recently forming. Numerous 
examples show that they are typically irregularly shaped and 
occasionally undergoing violent merger events. Cluster-sized 
dark-matter halos in simulations can often be well descr ibed as 
triaxial ellipsoids, but not as spheres ll.ling & ■Sutol2nn2h . 


At the same time, observations of galaxy clusters 
are often interpreted based on spherically-symmetric 
models in hydrostatic equilibri um. The beta model 
dCavaliere & Fusco-Femian^ Il976l) is still routinely be¬ 
ing used for analyses of the X-ray emission and also of the 
amplitude of the thermal Sunyaev-Zel’dovich effect. Given 
the rapidly improving quality and diversity of cluster data, it 
appears timely to search for an algorithm which avoids the 
assumption of spherical symmetry and allows the joint analysis 
of different types of cluster data. 


Several such algorithms have been proposed. IZaroubi et alJ 
suggested to base the reconstruction of axisymmet- 


ric, three-dimensional gravitational cluster potentials on the 
Fourier slice theorem, extrapolating Fourier modes into the 
“cone of ignorance”. They applied their techn ique to simu¬ 
lated d ata and sh o wed t hat it performs well dZaroubi et al.l 
I2OOII). iDore et alJ ll200lh followed a perturbative approach, 
and luee & Sutd ( 12004 ) proposed to adapt parameters of triax¬ 
ial halo models, all by combining different data sets such as 
X-ray, (thermal) Sunyaev-Zel’dovich (SZ) and gravitational- 
lensing maps. A simil ar method was applied to data by 
IPe Filippis et alJ ll200-5l) . An alternative approach based on 
the iterative Richardson-Lucy deconvolution was suggested by 


iReblinsk^ (l2000l) and iReblinskv & Bartelmannl(l200ll) . It aims 
at the gravitational potential, assumes only axial symmetry of 
the main cluster body, avoids extrapolations in Fourier space, 
and can easily be extended to include additional data sets. 

In this paper, we develop the latter algorithm further. 
However, aiming at the potential would require us to assume a 
relation between the gas distribution and the gravitational field, 
which would be most conveniently given by hydrostatic equi¬ 
librium. But even ignoring this common equilibrium assump¬ 
tion, it should be possible to reconstruct the three-dimensional 
distributions of intra-cluster gas density and temperature by a 
joint analysis of X-ray and thermal SZ data. 

We demonstrate here that this is indeed possible under 
the one simplifying assumption that the underlying three- 
dimensional distributions be axially (not spherically!) sym¬ 
metric. The inclination of the symmetry axis can be arbitrary. 
We introduce the algorithm in Sect. 2 and apply it to the ide¬ 
alised case of an analytically modeled, ellipsoidal cluster with¬ 
out substructure in Sect. 3. Results obtained first without, then 
with observational noise are highly promising; both the three- 
dimensional density and temperature distributions are accu¬ 
rately reproduced. Noise suggests smoothing, and we describe 
a suitable smoothing algorithm. 

We study the less-ideal case of a numerically-simulated 
galaxy cluster in Sect. 4. Here, axial symmetry is typically vi¬ 
olated by the main cluster body, and substructures are present 
which further perturb the symmetry. Yet, faithful reconstruc¬ 
tions are possible even in presence of realistic noise. 

Section 5 finally describes how inclination angles can be 
constrained using temperature maps, and Sect. 6 summarises 
and discusses our results. 
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2. The deprojection algorithm 


2.1. Deprojection of axisymmetric quantities using 
Richardson-Lucy deconvolution 

As iBinne y et alJ ( Il990l) poi nted out, Richardson-Lucy decon- 
volution iIlucvI1i97411i 994|) can be used to reconstruct an in¬ 
clined axisymmetric three-dimensional distribution of some 
physical quantity (j) from a two-dimensional map y/ of its pro¬ 
jection along the line-of-sight. In astrophysical applications, 
\j/ will be data obtained from observations, for example the 
X-ray flux, the lensing potential, or the Sunyaev-Zel’dovich 
decrement of an approximately axisymmetric galaxy cluster. 
Because of the assumed symmetry, (j) can be written as a func¬ 
tion of only two cylindrical coordinates R and Z, where we 
choose the symmetry axis as the Z-axis (see Fig. EJ- Then, R 
is the distance from the symmetry axis. The projection along 
the line-of-sight can be understood as a convolution of (j) {R,Z) 
with a kernel function P{x,y\R,Z), 

y/{x,y) <Pix,y,z) = JTtdR^j dZ (j){R,Z) P{x,y\R,Z). 

( 1 ) 

The kernel function for a given pair {R,Z) is non-zero only on 
the ellipse obtained by projecting the ring onto the sky which 
is defined by R and Z (see Fig.Q. It is derived in Appendix A 
of fRinnev et al .1 (fl 99f)l) and reads 


P{x,y\R,Z) 


5[(^-ZtanO^-(R^-^^)] 
n cos i 


( 2 ) 


where i is the inclination angle of the symmetry axis, defined 
as the angle between the symmetry axis and the line-of-sight; 
see Fig. [I] The kernel P{x,y\R,Z) satisfies the normalisation 
condition ^ ^ 

J dx J dy P{x,y\R,Z) = 1. (3) 

Assuming that the orientation of the symmetry axis is known 
and that one has a two-dimensional map of the projection y/, 
one can reverse the convolution us ing the iterative Ri chardson- 
Lucy deconvolution technique (see lLuc vll 97411 994 and solve 
for ^ as a function of R and Z. Starting with an initial guess 0o 
for (j) and using the Richardson-Lucy iteration scheme. 


(j)„+i(R,Z) = (j)„(R,Z) / dx Hdy Pix,y\R,Z), 

J-oo J-oo y/„{x,y) 


(4) 


one can obtain approximations of (j) with increasing quality. 
Here, y/„ is the projection along the line-of-sight of the approx¬ 
imation (j)„. If we plug Eq. 0 into Eq. 0, perform the y inte¬ 
gration, and use th e new coordinate a defined hy x = R cos a, 
we obtain (see also lBinnev et al .ll99(ilReblinskvl2000ll 


(j)„+i(R,Z) R^da y/(Rcosa^Zsini + Rsinacosi) 
^„{R,Z) Jo 2k y/,, (Rcosa,Zsin/-f-Rsinacosi) ’ 

where the integration follows the ellipse shown in Eig.^ 

Eor numerical reconstructions of axisymmetric three- 
dimensional distributions, we replace the integral in Eq. 0by 
a sum over points, which are distributed along the ellipse and 
equally spaced in a. To evaluate the sum, we need to find the 


P(x,y \R^Z) is non-zero here 




Fig. 1. Projection of an axisymmetric distribution of a physi¬ 
cal quantity. The ellipse at the top marks the region where the 
kernel function corresponding to the projection along the line- 
of-sight is non-zero for fixed R and Z. 


ratio yf/yfn at these points. We use two grids for the iteration, 
one in x,y-space for y/ and y/„, and one in R,Z-space for 
Eirst, we project <j)„ along the line-of-sight on the grid in x,y- 
space to find y/„. We do not use the kernel function for that, but 
perform a direct summation using a discretised version of the 
first equality in Eq. o . The projection integral is approximated 
by a sum over equally spaced points that cover a section of 
length Lj of the line-of-sight. This section is centred on the z- 
coordinate of the halo. Then, y/„ is obtained by 

Wn{xj,yk) = ■^Y^<f)n{R{xj,yk,zi),Z{xj,yk,zi)), (6) 

where Xj and yk are the x and y coordinates of the grid point 
{j,k). The zt are the z coordinates of the points used for 
the projection. The function ^„{R{xj,yk,zi),Z{xj,yic,zi)) is ap¬ 
proximated by bilinear interpolation from the values of (/)„ at 
neighbouring grid points in R,Z-space. Since we know y/ and 
have calculated the projection (//„ of 0„, we can find the ratio 
y//y/n at points on the ellipse by bilinear interpolation from 
neighbouring points of the x,y-space grid. This allows us to 
approximate the integral in Eq. (0 by a sum over Na points, 

(j)„+i{R,Z) 1 S (Rcosa,„,Zsin/-t-Rsina,„cos!) 

0„(R,Z) Wn (Rcosam,Zsin/ + RsinamCos/)’ 

(7) 

and find which completes the iteration step. Here a„, = 
2Km/Na- 
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2.2. Boundary efFects, artifacts and regularisation Then we set 


There is, however, a problem. Assume that corresponds to 
the height of the box shown in Fig.Q] and that the area covered 
by our map of xj/ corresponds to its top surface. To calculate 
V/n+i there, we have to know everywhere in the box. But 
for finding close to the corners of the box, we have to 
evaluate Eq. Q along ellipses that do not fit into the top sur¬ 
face of the box. This means that some of the points we have 
to sum over l i e outs ide our map of \j/ and v/„. As suggested 
by iReblinskvl ll2000ll . we replace xjf/xyn for those points by its 
value at the closest point at the perimeter of the map. This leads 
to some artifacts in the reconstruction of <f) for large R and Z, 
but yields very good results in the central region, which we are 
most interested in. 

To start the iteration, we have to choose a guess or prior ^o- 
We adopt the simplest choice of a flat or constant prior. We set 
its value so as to reproduce the average value (xj/) of the map 
xy, namely 00 = (r)/T^. 

The algorithm described above can be used to recon¬ 
struct axisymmetric three-dimensional distributions from two- 
dimensional maps of its projection along the line-of-sight. 
However, it runs into problems for strongly peaked distribu¬ 
tions such as the X-ray emissivity of a galaxy cluster. In order 
to illustrate that, we reconstructed the X-ray emissivity from 
an X-ray flux map, which we obtained by projecting the emis¬ 
sivity of an analytically modeled, axisymmetric cluster halo. 
The halo model is discussed in Sect. 0 For the projection, we 
chose an inclination angle of i = 70°. We performed the re¬ 
construction with a rather high number of n = 30 iterations. In 
the left panel of Fig. |2 we show the ratio between the recon¬ 
structed and the original X-ray emissivity. One can clearly see 
spike-shaped artifacts of the reconstruction. The angle between 
these spikes and the symmetry axis is equal to the inclination 
i. This means that the ellipses corresponding to R and Z values 
of points in the spikes pass directly through the halo centre in 
the map of xj/. 

Richardson-Fucy deconvolution reproduces large struc¬ 
tures quickly, while it converges slowly to small scale struc¬ 
tures such as the peak at the halo centre (see IFucvI I1974 
I 1994 I) . This means that, when starting with a flat prior, xy/xj/n 
can be quite large close to the centre even after several it¬ 
erations. Thus, when we evaluate 0 for points further out 
whose ellipses pass through the halo centre, we And ratios of 
0„+i(R,Z)/0„(R,Z) which are too high, and the spike-shaped 
artifacts form. They appear already after the first few iterations 
and are very stable. In the left panel of Fig. |2l we show them 
after 30 iterations, and it would take several hundred more iter¬ 
ations until they slowly disappear. 

To prevent the formation of such artifacts, we use a regular¬ 
isation scheme. First, we calculate an average {xj//xj/„) for the 
points used in the sum in Eq. 0 , which is defined by 

{V/Vn)iR,Z} = 

1 ^ . / iy(Rcosam,Zsin/ + RsinamCosi) 

= —) mm —7 - 7)10 • 

Vy/„(Rcosam,Zsin/-t-Rsina,„cosij / 

(8) 


Cn{R,Z) 


1 f cos Om , Z sin / -f R sin cos i) 

— > mm I —---, 

Na ^ xj/n(R cos Urn,Zsini + R sm a,,, cos i) 

l.2{xy/xyn){R,Z)), (9) 


and use 

=max(c„(R,Z), 0.25 ), (10) 

to calculate 0„+i(R,Z). This regularisation of the iteration 
scheme supresses the formation of the spike-shaped artifacts. It 
limits the impact of sharp peaks in xj/ on points that are far away 
from the corresponding peaks in 0 by using an upper limit for 
xy/xj/n in Eqs. (|^ and (|3. The effectiveness of the regularisa¬ 
tion scheme is not very sensitive to the exact numerical values 
of the upper limits, which we chose by trial and error, as long 
as they suppress sharp peaks and are not too restrictive to allow 
convergence in a reasonable number of iterations. The lower 
limit for the correction factor in Eq. ( I10> is introduced just to 
make sure that 0„+i does not change its sign or become very 
small in the first few iteration steps, which would potentially 
cause problems later. 

We repeated the reconstruction of the X-ray emissivity us¬ 
ing this regularisation. In the right panel of Fig. |2l we show 
the ratio of the reconstructed to the original emissivity after 
30 iterations. The spikes that are present in the left panel have 
almost disappeared. The ratio is close to unity everywhere in 
the region shown, except very near the halo centre where grid 
resolution and the slow convergence to small-scale structures 
becomes a problem. Apart from that, the deprojection works 
very well. The errors are usually smaller than 1%. 


2.3. Inclination angle 

So far, we have assumed that the orientation of the symmetry 
axis is known. In reality, when applying this algorithm to ob¬ 
servations, this will not be the case. However, the orientation 
of the symmetry axis in the plane of the sky can directly be ob¬ 
tained from the map xy. Its inclination i on the other hand can 
be found in the following way. First, we repeat the reconstruc¬ 
tion as described above, assuming several different plausible 
values of i and using a Axed number Ait of iterations. Then, 
we compare the maps xyj^.^^, corresponding to the reconstructed 
distributions, to the original (or observed) map xy and And the 
value of i for which it fits best. This can for example be done 
by minimising the penalty function 

ivixj,yk) - xyN^^{xj,yk)f 

where, depending on the shape of the distribution one wants to 
reconstruct, it may be favourable to sum only over points in the 
central region of the map. 
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without legularisation with regulaiisation 



0.8 1 1.2 


Fig. 2. Ratio of the reconstructed to the original X-ray emis- 
sivity after n — 30 iteration steps, The halo centre is at 

the centre left of each plot. Panel (a) shows the ratio obtained 
without regularisation. One can clearly see the spike-shaped 
artifacts of the reconstruction. Panel (b) shows the ratio for the 
reconstruction including the regularisation. It is close to unity 
everywhere except very close to the centre, where the algorithm 
converges slowly and grid resolution becomes a problem. 


2.4. Reconstruction of density and temperature of the 
ICM from combined X-ray and thermal 
Sunyaev-Zel'dovich effect observations 


where the integration also follows the line-of-sight, Icb is 
Boltzmann’s constant, nig the electron mass, c the speed of 
light, Gt the Thomson cross section, and rig the electron num¬ 
ber density at the position passed by the light ray at time t. For 
fully ionised gas with constant metallicity, this is also propor¬ 
tional to \j/sz defined by 

\l/sz = jdzpT. (14) 

V^x-ray ^nd xj/sz Can be obtained from observations. For recon¬ 
structing the ICM temperature and density, we start from some 
initial guess po{R,Z) and Tq{R,Z). In analogy to Eq. we 
use discrete approximations of Eqs. O and 113, 


Lj 


VX-my,n{xj,yk) = J^Y.Pn{^^^j^yk,Zl),Z{xj,yk,Zl)) 


z/=l 


X X{T{R{xj,yk,Zi),Z{xyyk,Zi)),S’), (15) 


Lj 


Nz 


WszA^^yk) = -^Y^Pn{R{xj,yk,zi),Z{xj,yk,zi)) 
X T{R{xj,yk,Zi),Z{xj,yk,Zi)), 


(16) 


to obtain i/x-ray,o Wsz.o- In analogy to Eq. Q, we define 
Zr. 1 ^ V6c-ray COS am,Z sin i-fR sin cos/ 

CX-ray,«(f^,Z) =— 2^ 


Na i/x-ray,n {R COS , Z sin / +R sin am cos /) ’ 

(17) 

^ V V^sz(f^cosa,„,Zsin/-t-RsinamCos/) 

J.J / . / ’-y ' ‘ \ n ' ‘X ’ 

^1 rsz, n cos a™, Z sm! 4-R sm a,„ cos!) 

(18) 


So far, we have discussed how to reconstruct a single three- 
dimensional distribution of a physical quantity from a sin¬ 
gle two-dimensional map of its projection along the line- 
of-sight. However, one can obtain additional information 
by combining maps from different observati ons (see also 
Reblinskv & BartelmanrJ2Q0lUReblinskvi2003l . Here, we re¬ 
construct several physical quantities at the same time by com¬ 
bining different observations that depend on these quantities. 
Specifically, we will show how to obtain three-dimensional dis¬ 
tributions of the density and temperature of the intra-cluster 
medium (ICM) in axisymmetric cluster halos by combining 
observations of the (thermal) X-ray emission and the (ther¬ 
mal) Sunyaev-Zel’dovich effect. The X-ray surface brightness 
is taken to be proportional to 

XI/X-rzzy = Jdzp^XiT,^), (12) 

where p and T are the gas density and temperature, respec¬ 
tively. The integral extends along the line-of-sight. The cooling 
function X{T,3f) depends on the temperature and the metal¬ 
licity and is assumed to be constant here. The (thermal) 
Sunyaev-Zel’dovich decrement or increment is proportional to 
the Compton y-parameter 


and use the iteration scheme 


plX[Tn,3Z) 


Pn+\Rn-\-\ 

PnPn 


^X-ray, 


CSZ,«- 


(19) 

( 20 ) 


In order to find the next iterative approximation for the den¬ 
sity and the temperature, we have to solve equations J19t and 
(El for p„+i and T„^i. To include line emission in the cooling 
function, one could tabul ate X{T„. .^), e.g. using the software 
package XSPEC (see e.g. lArnaudll996l) for a specific emission 
model and metallicity, and solve the equations above numeri¬ 
cally. However, since the main focus of this paper is to demon¬ 
strate the deprojection algorithm, we use the simple assumption 
of continuous bremsstrahlung. 


A(7’, 2f) = Vf. 

Then, we obtain from ( I19> and EqJ, 

2/3 

r 

Pn+l = 


‘'X-ray, n 


( 21 ) 


X-ray,« 

^1/3 

Pttj 

(22) 

‘'SZ,« 



4/3 



^SZ,n 

,2/3 

■Tn. 

(23) 


y= / - lyneOTCdt, 

J nigC^ 


(13) 


Tn+l = 
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Note that for evaluating Eq. (O and Cl, we also use the reg- 
ularisation introduced in Sect. 12.2| and expressed by Eqs. (H- 
d- Again, for points that lie outside the map of i/x-ray and 
Xj/sz, the ratios V^x-ray/VA:-ray,« and xi/sz/Wsz,n are approxi¬ 
mated by their values at the closest point at the perimeter of the 
map. In the next two sections, we shall apply this deprojection 
algorithm to axisymmetric analytic halos and to numerically 
simulated cluster halos, and discuss its performance. 


3. Deprojection of analytic halos 

The analytic halo model which we use to test the deprojection 
algorithm is based on the NEW density profile. However, for 
the deprojection to be non-trivial, we prefer to have ellipsoidal 
iso-density surfaces. We therefore assume that the density of 
the ICM is a function of 


r = 





(24) 


where is a scaling radius perpendicular to the symmetry 
axis, and Z* is a scaling distance along the axis. The density 
of the hot cluster gas is then taken to be 


P = 


Ps 

(er + r)(l +r)2’ 


(25) 


which differs from an NEW halo not only by its ellipsoidal 
shape but also by the small constant Sr = 0.001 introduced to 
ensure that the density does not diverge for r ^ 0. This diver¬ 
gence would cause problems in numerical calculations. Eor the 
gas temperature we use a phenomenological description that 
roughly corresponds to the temperature profiles found in the 
simulated cluster sample that we will use in Sect.^] Namely, 
we set 

r = (26) 


where /(r) = tanh(3(r— 1)) is —1 for r ^ 1 and -f 1 for r ^ 1. 
The values of the parameters = 3 x lO^/j^*M0/(/j^^kpc)2, 
Zs = 500/i^'kpc, Rs = 300/i^'kpc, and Tmax = 12keV, which 
we use correspond to the gas component of a massive galaxy 
cluster. Here, h is the reduced Hubble constant, which we set 
to 0.7. 


3.1. Deprojection of analytic halos without 
observational noise 

Having chosen an inclination angle i, we project the analytic 
halo described above on a 128 x 128 grid with a sidelength of 
1.5/i^^Mpc and obtain the X-ray and Sunyaev-Zel’dovich ef¬ 
fect maps, i/x-ray and y/sz- We use the algorithm discussed in 
Sect. 12.41 with these maps to reconstruct the gas density and 
temperature. In Eig.|3 we compare the results of the deprojec¬ 
tion with the original density and temperature of the analytic 
halo for an inclination angle i = 70° and after n = 20 iterations. 
The inclination was assumed to be known in performing the 
deprojection. The left and right panels show the density and 
temperature ratios, p,j/p and T„/T, repsectively. The star-like 
pattern of the plots maps the ranges of the R and Z coordinates 
occurring in the simulation box used for the reconstruction (see 


Fig.lD. In the central region of the cluster, the reconstruction 
works very well. Errors are of the order of 1%. Despite the reg- 
ularisation, one can still see some remains of the spike-shaped 
artifacts discussed in Sect. 12.21 Eor large R or Z values, close 
to the star-shaped boundary of the plots, the quality of the re¬ 
construction decreases. This is not at all surprising because the 
ellipses along which Eqs. o and (Hill must be evaluated to 
reconstruct density and temperature at those points lie mostly 
outside of the maps of xj/x-ay and xf/^z- 


density ratio temperature ratio 



0.95 1 1.05 


Fig. 3. Ratio of the reconstructed to the original density and 
temperature for the analytic halo after n = 20 iterations. An in¬ 
clination angle of i = 70° was chosen and assumed to be known 
in performing the reconstruction. Panel (a) shows the density 
ratio Pn/p and panel (b) the temperature ratio T„/T. In the cen¬ 
tral region, the errors are of the order of 1%. We plot the ratios 
for all R and Z values possible within the box used for the re¬ 
construction (see Eig.Q. This causes the star-like shape of the 
perimeter of the plot in the R, Z plane. Close to that boundary, 
at large R or Z values, the ratios can significantly differ from 
unity. This is, however, expected because the ellipses used in 
the reconstruction of p and T at those points lie mostly outside 
the maps of i/x-ray and y/sz- 


Note that the quality of the reconstruction also depends on 
the inclination of the halo’s symmetry axis. Of course, best re¬ 
sults are achieved when the symmetry axis is perpendicular to 
the line-of-sight. Then the assumption of axial symmetry con¬ 
tains the most information. If, on the other hand, the symmetry 
axis is parallel to the line-of-sight, the axial symmetry just cor¬ 
responds to the circular symmetry of the maps i/x-ray and xysz 
and does not yield any useful additional information. Eigure|4] 
illustrates this inclination dependence. It shows the volume- 
weighted root mean square (RMS) relative errors of the recon¬ 
structed gas density and temperature, computed within a sphere 
of radius 500/i^'kpc around the halo centre. Again, the knowl¬ 
edge of the inclination angle i was used in the deprojection. 
An accuracy of 1% or better is achieved for about two thirds 













6 


Puchwein and Bartelmann: Three-dimensional reconstruction of the intra-cluster medium 


of the analytic halos in a randomly oriented sample. However, 
halos that happen to have a very small inclination angle i are 
necessarily poorly reconstructed. 

inclination i [°] 


10 20 30 40 50 60 70 80 90 



Fig. 4. Dependence of the quality of the deprojection of ana¬ 
lytic halos on the inclination of the symmetry axis. We show 
the volume-weighted RMS relative errors {pn — p)/p and 
{Tn — T)/T as functions of the inclination angle i after n — 20 
iterations. The averages were computed within a sphere of ra¬ 
dius 500/i^'kpc around the halo centre. The quantity 1 — cos(/) 
shown on the abscissae is chosen such as to have a flat number- 
density distribution for randomly oriented halos. The inclina¬ 
tion angle i was assumed to be known in performing the depro¬ 
jection. 


3.2. Deprojection of analytic halos including 
observational noise 

So far we have not considered noise that will be present in any 
real X-ray or Sunyaev-Zel’dovich effect observation. We will 
now discuss the impact it has on the reconstruction of ICM 
densities and temperatures. 

We model the noise in X-ray observations as follows. First, 
we calculate for each pixel {j, k) of the halo’s X-ray map i/x-ray 
the number of photons {Nyj^j) expected from bremsstrahlung, 

which is proportional to {Nyj^k) "yT’ where the 

sum extends along the line-of-sight represented by the pixel 
( 7 ,k), and Ei is the exponential integral function. Frnin is a 
lower energy cutoff which is necessary because the number 
of photons emitted is infrared divergent. We choose Emin = 
0.23keV, which is a reasonable lower limit for the photons 
from galaxy clusters observed in current X-ray experiments. 
Next, we normalise the numbers of expected photons such that 
they sum up to Y.i.k{f^yi,k) = 10"^ on the entire map. For each 
pixel (y'jk), we then set the actual number of photon counts 
Nyj f^ to a value drawn from a Poisson distribution with expec¬ 
tation value {Nyj_k). We then add the noise to the map i/x-ray 
by multiplying V^x-iay;.t: with Nyj.k/{N'Yj,k) for pixels. 

For the Sunyaev-Zel’dovich effect, we add noise cor- 
responding to future ALMA Band 3 observations (see 
iButler & Woottenlll999l) . In Band 3 (84 to 116 GHz) and in 
its compact configuration, ALMA will be able to achieve a 


temperature sensitivity of 50/rK at a spatial resolution of ^ 3 
arcsec in about four hours of observation. At an assumed halo 
redshift of 0.3, this resolution corresponds to the angular size 
chosen for the pixels of the map i/sz- We convert the temper¬ 
ature sensitivity cited above to an error of ¥sz- Then, for 
each pixel, we add noise obtained from a normal distribution 
with standard deviation to y/sz- 

Richardson-Lucy deconvolution has the nice property of 
approximating large-scale features quickly and small-scale 
noise slowly. Yet, it turns out that smoothing the noisy maps 
VN-ray and \gsz before using them in the deprojection im¬ 
proves the results considerably. We use the following smooth¬ 
ing scheme. Lor the X-ray observations, we assume that in ad¬ 
dition to the map Xj/x-ay we also know the photon counts Nyj ^ 
for all pixels. We then calculate for each pixel {j,k) a radius 
hj i^ so that we have a fixed number of 100 photons inside a cir¬ 
cle with radius hj^ around that pixel. After that we redistribute 
the value Yx-ray j.k of each pixel on the grid with a smoothing 
kernel of width hj j centred on that pixel. This greatly reduces 
the fluctuations in the map V6c-ray caused by photon noise. In 
the remainder of the paper we will call this first step of the 
smoothing scheme “photon-noise smoothing”. 

Lor the smoothing kernel, we take the line-of-sight projec¬ 
tion of the cubic spline SPH kernel W {r,h) defined in Appendix 
A of ISpringel et alJ (120011) . It is well suited for this purpose 
and allows us to use the same routine for smoothing here and 
for projecting the numerical SPH halos we use in Sect. 0] Lor 
axisymmetric halos, the projection should be symmetric about 
the projected axis. However, the symmetry is broken here by 
noise. We restore it before performing the deprojection. Since 
we have oriented the grid such that it is parallel to and cen¬ 
tred on the projected symmetry axis, we can do that by replac¬ 
ing V6c-rayj,t: xgx-myj.k arithmetic mean. Here, 

Agrid = 128 is the dimension of the grid. We symmetrise i/sz in 
the same way. After that we perform one more smoothing op¬ 
eration on i/x-ray and i/sz to further reduce fluctuations caused 
by noise. 

In numerically simulated halos, which we will discuss later, 
this will also suppress the effect of subclumps. Since we do 
not want to smooth out the peaks in the halo core, we choose a 
smoothing length h that depends on the distance r from the halo 
(or map) centre, namely /j =/imax(l - W(r,rmax)/W(0,rn,ax))■ 
It is zero in the centre of the map and continually increases to 
h — 375/;^ *kpc at a radius equal to r^ax = 150h^ *kpc or larger. 
This yields smallest RMS errors. Once we have calculated h 
for each pixel, we smooth the maps of V6t-ray and VAz with 
the projection of the SPH smoothing kernel mentioned above 
and with the position-dependent smoothing length h. Note that 
roughly 80% of a pixel’s value is redistributed within a circle 
of radius hjl. We refer to this second step of the smoothing 
scheme as “radius-dependent smoothing”. 

After degrading the analytic halo with noise and applying 
the smoothing scheme described above, we perform the itera¬ 
tive deprojection. The results after n = 5 iterations are shown 
in Pig. 13 An inclination of i = 70° was chosen and assumed to 
be known in the deprojection. Again, the left panel shows the 
ratio of the reconstructed to the original density, and the right 
panel the corresponding temperature ratio. Average errors in 
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the central region are of the order of 5% to 10%. As expected, 
further outside, where the signal-to-noise ratio becomes small 
and the ellipses used for the reconstruction lie mostly outside 
the maps of iffx-my and if/sz, the errors are substantially larger. 
Note that at locations where we obtain a too low density, we 
usually hnd a too high temperature and vice versa. This hap¬ 
pens because the algorithm minimises the deviations of the re¬ 
constructed from the original X-ray and SZ-effect maps. 



Fig. 5. Ratio of the reconstructed to the original density and 
temperature for the analytic halo with observational noise after 
n = 5 iterations. An inclination of i = 70° was chosen and as¬ 
sumed to be known in performing the reconstruction. “Photon- 
noise” and “radius-dependent” smoothing were applied. Panel 
(a) shows the density ratio Pn/p, and panel (b) the temperature 
ratio Tn/T. In the central region, the errors are of the order of 
5% to 10%. 


In Fig. |6l we show density and temperature profiles of the 
original analytic halo, of the halo reconstructed from maps 
without observational noise, and of the halo reconstructed from 
smoothed maps which contain observational noise. The recon¬ 
structed halos are the same as shown in Figs.|3]and|3 Without 
noise, both the temperature and the density prohles are repro¬ 
duced very well. With noise, we can still reproduce density 
prohles with an accuracy of a few percent. The errors in the 
temperature prohle are somewhat larger. Deviations are mainly 
caused by the noise, but some are also artifacts of the smooth¬ 
ing scheme we apply. Especially the too high temperature near 
~ 200kpc and the too low temperature near ~ 400kpc 
are a consequence of “radius-dependent smoothing”. On the 
other hand, without such smoothing we would approximately 
double the etTors in the density and temperature reconstruc¬ 
tions. 



Fig. 6. Density and temperature prohles of the original and the 
reconstructed analytic halos. The upper panel shows the density 
(falling curves, left axis) and the temperature prohles (rising 
curves, right axis) of the original analytic halo, the halo recon¬ 
structed without observational noise (and without any smooth¬ 
ing), and the halo reconstructed from maps with observational 
noise to which the complete smoothing scheme was applied. 
The lower panel shows the prohle of the ratio of the recon¬ 
structed density p„ to the original density p. The number of 
iterations used was n = 20 in the case without noise and n = 5 
in the case with noise. An inclination of i = 70° was chosen 
and assumed to be known in the reconstruction. 


4. Deprojection of numerically simulated cluster 
halos 

So far, we have demonstrated the performance of the algorithm 
with axisymmetric analytic halos. We were able to reconstruct 
their three-dimensional density and temperature distributions 
from synthetic X-ray and (thermal) Sunyaev-Zel’dovich effect 
observations. Real galaxy clusters, however, are hardly per¬ 
fectly axisymmetric. We will study in this section whether they 
nonetheless allow accurate density and temperature reconstruc¬ 
tions with the deprojection algorithm proposed in Sect. 12.41 We 
use a sample of four numerically simulated galaxy clusters to 
investigate into this question. 

The simu lations were ca rried out by Klaus Dolag with 
GADGET-2 ISmngd E^, a new v ersion of the paralle l 
TreeSPH simulation code GADGET dSoringel et alJliooIl) . 
The cluster regions were extracted from a dissipation-less (dark 
matter only) simulation with a box size of 479/1^^ Mpc of 
a flat ACDM mode l with = 0.3, h = 0.7, ds = 0.9 (see 
IVoshida et al ■I2OOII) . They were re-simulated with higher reso- 
lution using th e “Zoomed Initial Conditions” (ZIC) technique 
(lTormerill997h . Gas was introduced into the high-resolution 
region by splitting each parent particle into a gas and a dark 
matter particle, which were then displaced by half the mean 
inter-particle distance, such that the centre-of-mass and the mo¬ 
mentum were conserved. The mass ratio of gas to dark matter 
particles was set to obtain Clh — 0.04. The hnal mass resolution 
was otdm = 1-13 X IO^/j^'Mq and mgas = 1.7 x IO^/i^'Mq 
for dark-matter and gas particles within the high-resolution re¬ 
gion, respectively. The simulations we use follow the dynamics 
of the dark matter and the adiabatic evolution of the cluster gas. 
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but they ignore radiative coolin g. They are described in more 
detail in IPuchwein et al.l (l2005l) and lOolag et alJ (l2005l) . 

Our deprojection algorithm requires a symmetry axis, 
which real clusters do not generally have. We thus need to 
choose an axis around which the numerical clusters have at 
least a high degree of symmetry. We do this by calculating 
the inertial tensor of the cluster gas inside a sphere of radius 
500/;^* kpc around the cluster centre and finding its eigenvec¬ 
tors vi, V 2 , V 3 and eigenvalues ej >e 2 > £ 3 . We choose the sym¬ 
metry axis through the cluster centre and parallel to the eigen¬ 
vector V 3 with the smallest eigenvalue if ei/e 2 < ^ 2 /^ 3 ^ or par¬ 
allel to the eigenvector vi with the largest eigenvalue otherwise. 
This means that, if two eigenvalues are very similar, we choose 
the axis parallel to the eigenvector corresponding to the third 
eigenvalue. 

4.1. Deprojection of numerical halos without 
observational noise 

Having chosen a fiducial “symmetry” axis and a line-of- 
sight, we can produce synthetic maps of X-ray and Sunyaev- 
Zel’dovich effect observations. For that purpose, we use the 
simulated clusters at a redshift of z = 0.3 and project them 
along the line-of-sight. At z = 0.3, the cluster sample spans 
a mass range between 8 x lO'"^ and 1.8 x lO'^h^^M©. 

For now, we do not add any observational noise to the 
maps. However, the clusters contain substructures which break 
axial symmetry and lead to artifacts in the density and tem¬ 
perature reconstructions. Thus, depending on the amount of 
substructure present in a cluster, it may still be favourable to 
use “radius-dependent smoothing” on the X-ray and Sunyaev- 
Zel’dovich effect maps prior to reconstruction. In Figs. □andE] 
we show the results of the deprojection without any smoothing, 
and using “radius-dependent smoothing”. 

The density reconstruction in the central region reaches an 
accuracy of about 10% in both cases. For the temperature re¬ 
construction and the density reconstruction at large r, we obtain 
somewhat better results without smoothing for this rather sym¬ 
metric cluster. Note, however, the hyperbolically shaped arti¬ 
facts in Fig. 0 which are produced by substructure clumps in 
absence of smoothing. They appear at those R and Z values 
which correspond to the line-of-sight passing through such a 
clump. The spike-shaped artifacts discussed in Sect. l2.2l were a 
special case of the artifacts found here. For most of the hyper¬ 
bolae in the left panel of Fig. 0 one can also see the position 
of the clump that produced it in darker colours. The hyperbolae 
pass right through them. 

As one can see in Fig.| 8 ] “radius-dependent smoothing” re¬ 
moves the hyperbolic artifacts. The subclumps, however, still 
appear in darker colours in the density ratio map, which means 
that the reconstructed density there is too low. However, this 
is entirely expected and inevitable, because they violate axial 
symmetry and thus cannot be faithfully reconstructed with this 
deprojection technique. By smoothing, we essentially remove 
the subclumps from the maps and reconstruct the density and 
temperature of the main halo without them. 


Unfortunately, “radius-dependent smoothing” also affects 
the density and temperature profiles. This can be seen in the 
“rings” around the halo centre in the right panel of Fig. |8] It 
is further illustrated in Fig. which shows the density and 
temperature profiles of the original cluster g51, after deprojec¬ 
tion without noise but with “radius-dependent smoothing”, and 
after deprojection without noise and without smoothing. For 
r > 300 / 1 ^ *kpc, the reconstruction without smoothing yields 
more accurate density and temperature profiles. In addition, the 
profiles for deprojections from maps including observational 
noise are shown. They will be discussed in the next section. 

Reconstructions along different lines-of-sight and of the 
three other clusters in the sample gave similar results. For the 
most asymmetric halo, the errors were larger by factors of 1.5 
to 2 compared to the reconstruction of g51 presented above. 
We can thus conclude that, although clusters are not strictly 
axisymmetric and contain substructure, it is possible to apply 
the deprojection method proposed in Sect. l2.4l and successfully 
reconstruct three-dimensional density and temperature distri¬ 
butions of the cluster gas. 


density ratio temperature ratio 



0,8 1 1,2 


Fig. 7. Reconstruction of the simulated cluster g51 without 
noise and smoothing. The ratios of the reconstructed to the 
original density and temperature are shown. The deprojection 
was done with n — 5 iterations. An inclination angle of i = 68° 
between the line-of-sight and the principal inertial axis of the 
cluster gas was chosen and assumed to be known in the recon¬ 
struction. Panel (a) shows the density ratio Pn/p and panel (b) 
the temperature ratio Tn/T. In the central region, the errors are 
of the order of 10%. 


4.2. Deprojection of numerical halos including 
observational noise 

In Sect. 13.21 we studied the impact of observational noise in 
the X-ray and Sunyaev-Zel’dovich effect maps on the quality 
of the density and temperature reconstruction. We will now do 






















Puchwein and Bartelmann: Three-dimensional reconstruction of the intra-cluster medium 


9 




500 


500 


-500 


-500 


R[h *kpc] 

(a) 


1000 

(b) 


density ratio 


1000 


temperature ratio 


1000 


R[h *kpc] 500 


density ratio 


temperature ratio 



Fig. 8. Reconstruction of the simulated cluster g51 without 
noise but with “radius-dependent smoothing”. The ratios of 
the reconstructed to the original density and temperature are 
shown. The deprojection was done with n = 5 iterations. An 
inclination angle of i = 68° between the line-of-sight and the 
principal inertial axis of the cluster gas was chosen and as¬ 
sumed to be known in the reconstruction. Panel (a) shows the 
density ratio p„/p and panel (b) the temperature ratio Tn/T. In 
the central region, the errors are of the order of 10%. 


Fig. 9. Reconstruction of the simulated cluster g51 with noise 
and and the complete smoothing scheme applied. The ratios 
of the reconstructed to the original density and temperature 
are shown. The deprojection was done with n = 5 iterations. 
An inclination angle of i ~ 68° between the line-of-sight and 
the principal inertial axis of the cluster gas was chosen and as¬ 
sumed to be known in the reconstruction. Panel (a) shows the 
density ratio p„/p and panel (b) the temperature ratio T„/T. In 
the central region, the errors are of the order of 15%. 


the same for the numerically simulated cluster halos using the 
same noise model, namely Poisson noise corresponding to 10^ 
observed photons for the X-ray maps and a noise level expected 
for a four-hour ALMA Band 3 observation for the Sunyaev- 
Zel’dovich effect maps. We also use the smoothing scheme de¬ 
scribed there. 

We show the results of the reconstruction in Fig.|3 Again, 
the left panel shows the ratio of the reconstructed to the origi¬ 
nal density, and the right panel the corresponding temperature 
ratio. In the central region we achieve an accuracy of about 
15%. Without “radius-dependent smoothing”, errors would be 
larger by roughly a factor of 1.5 or even more next to the 
halo centre. However, if one is mainly interested in density 
and temperature profiles it may still be favourable to leave the 
“radius-dependent smoothing” step away. Although the errors 
are larger without “radius-dependent smoothing”, they are less 
biased with respect to the distance from the halo centre and 
cancel better when averaging over spherical shells around it, 
especially at large radii. Figure^Jshows the profiles obtained 
with and without “radius-dependent smoothing”. 

We still need to discuss when the iteration used in the den¬ 
sity and temperature reconstructions should best be stopped. 
Figure illustrates the dependence of the quality of the re¬ 
construction on the number of iterations used. More precisely, 
it shows the relative volume-weighted RMS error of the density 
reconstruction within r = 500 * kpc as a function of the num¬ 
ber of iterations and for different deprojection schemes, namely 
for the deprojections of the analytic halo and the numerically 


simulated cluster g51 discussed above and shown in Figs.|3 
Elia and la after n — 20 or n — 5 iterations. The quality of the 
reconstruction improves quickly during the first roughly five it¬ 
erations (first ten for the analytic halo without noise) and then 
levels off. In addition, we show the quality of the reconstruc¬ 
tion of g51 from maps with noise but without using “radius- 
dependent smoothing”. In this case, small-scale noise in the 
maps is not sufficiently suppressed. The best reconstruction is 
found after five iterations. Then, the quality decreases again 
because the algorithm starts to approximate small-scale noise. 
Thus, unless a halo is very smooth and axisymmetric, such as 
the analytic halo without noise, we find that the quality of the 
reconstruction does not significantly increase after n = 5 iter¬ 
ations and may even decrease if small scale fluctuations due 
to noise are not efficiently suppressed. Thus, we conclude that 
it is favourable to use this number of iterations for the depro¬ 
jection of simulated and real galaxy clusters. Alternatively, one 
could control the reproduction of small-scale fluctuations with 
a formal regularisation scheme, such as provided by maximum- 
entropy methods. 

5. Finding inclination angles 

In all deprojections of analytic and numerical clusters pre¬ 
sented above, we have assumed that the inclination angle of the 
“symmetry” axis is known beforehand. This will usually not be 
the case for real observed clusters. In Sect. 12.31 we sketched 
how one can obtain inclination angles by comparing the maps 
obtained by projecting reconstructed halos to the original maps. 
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Fig. 10. Gas density and temperature profiles of the original 
and the reconstructed cluster g51. The upper panel shows the 
density profiles (falling curves, left axis) and the tempera¬ 
ture profiles (rising curves, right axis) of the original cluster, 
the cluster reconstructed without observational noise but with 
“radius-dependent smoothing”, reconstructed without observa¬ 
tional noise and without any smoothing, reconstructed from 
maps with observational noise and the complete smoothing 
scheme applied, and reconstructed from maps with observa¬ 
tional noise but without “radius-dependent smoothing”. The 
lower panel shows the profile of the ratio of the reconstructed 
density p„ to the original density p. The number of iterations 
used was n = 5 in all cases. An inclination angle of i « 68° was 
chosen and assumed to be known in the reconstruction. 



Fig. 11. Dependence of the quality of the density reconstruc¬ 
tion on the number of iterations used. We show the volume- 
weighted relative RMS error ^ ^ within spheres with r = 
500 kpc radius and for different deprojections of our an¬ 
alytic halo model and of the numerically simulated cluster g51. 
The figure legend explains which halo, whether or not noise, 
and what kind of smoothing were used for the different recon¬ 
structions shown. Note that “photon-noise smoothing” of the 
X-ray maps is always and only used for maps with noise. 


One would reconstruct a cluster assuming different values for i 
and find the value for which the maps match best. In principle, 
we could generalise this approach to the X-ray and Sunyaev- 
Zel’dovich effect maps used for the density and temperature 
deprojections. That is, we can compare the original maps V6c-ray 


and \j/sz to the maps V6c-iay,Vi, and i/sz. 7 Vit^ which correspond to 
the reconstructed halo after a fixed number Ait of iterations and 
for different values of the inclination angle i used for the recon¬ 
struction. 

We did this for the analytic halo model and for our sample 
of numerical clusters and used different inclination angles i' for 
projecting these halos to obtain the original maps. However, the 
minima in the penalty function (see Eq. GB) are not well de¬ 
fined. They are very broad and not always centred on i = 
Even for the analytic halo without observational noise, it is 
hardly possible to find the correct axis inclination in this way. 
As one can see from Eqs. o and (CD, the iterative correc¬ 
tions of the deprojection algorithm are determined from the 
deviations of the X-ray and Sunyaev-Zel’dovich effect maps, 
and the deviations are thereby minimised. Unfortunately, this 
still works remarkably well when choosing a wrong inclination 
angle i ^ i' for the deprojection. Thus the X-ray and Sunyaev- 
Zel’dovich effect maps are still reproduced well in this case, 
although the errors of the density and temperature reconstruc¬ 
tions increase significantly. 

We tried to limit the ability of the deprojection algorithm to 
reproduce maps well even when the inclination angle is wrong 
by reducing its degrees of freedom. Eor doing so, we used a 
variant of the algorithm that only reconstructs the density and 
uses a constant but adjustable temperature. This of course also 
limits the accuracy of the reconstruction for the correct incli¬ 
nation angle. Thus, the results of comparing the reconstructed 
X-ray and Sunyaev-Zel’dovich effect maps to the original ones 
for finding the inclination angle were not significantly better. 

On the other hand, leaving the deprojection algorithm as 
described in Sect. 12.41 but using additional independent infor¬ 
mation for finding the inclination angle of the halo, seems to 
be more promising. Eor the deprojection, we use maps of the 
X-ray flux of the clusters, but so far we do not use any spec¬ 
tral information from the X-ray observations. In Eigs. 1121 and 
[Tsi we assume that in addition to the X-ray flux maps we also 
have maps of the emission-weighted temperature Tew. We re¬ 
construct the analytic halo and the numerically simulated clus¬ 
ter g51 from X-ray flux and Sunyaev-Zel’dovich effect maps as 
above, but then compare the original emission-weighted tem¬ 
perature map to one we obtain by projecting the reconstructed 
halos. We repeated this for different inclination angles i', cho¬ 
sen for projecting the original maps, and i, chosen in the recon¬ 
struction. 

Eigure El shows the RMS relative error of the recon¬ 
structed emission-weighted temperature maps for the analytic 
halo without and with noise. Ait = 20 and no smoothing were 
used without noise, and Ait = 5 and the complete smoothing 
scheme were used with noise. The RMS error was computed 
within a radius of 500 kpc around the map centre and is 
shown for inclinations of i' = 40° and i' =70° of the original 
halo. As desired, the minima of the error curves are at the cor¬ 
rect locations i « 

Note that the curves are only shown for i between 0° and 
90° because Vfx-ray and i/sz and hence the whole deprojection 
algorithm is insensitive to what is the front and what is the back 
side of the cluster. We thus get the same reconstruction and 
the same errors for deprojections which adopt inclination an- 
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gles i and 180° — i. There is no way to distinguish these cases 
from the X-ray, thermal Sunyaev-Zel’dovich effect and temper¬ 
ature maps. The error curves are thus symmetric about i = 90°. 
Also note that for the halo with noise, we added observational 
noise only to the X-ray flux maps and Sunyaev-Zel’dovich ef¬ 
fect maps that were used for the reconstruction, but not to the 
Tew maps which we use for finding the inclination angle. We do 
not mimic observational noise in the temperature maps because 
it can only be realistically modeled when considering instru¬ 
ment response and line emission. In addition to the error curves 
of the temperature maps, we also show in Fig.^lthe volume- 
weighted, relative RMS errors of the density reconstructions in 
the central kpc. As expected, the reconstruction works 

best for i « i'. 

Figure^]shows similar quantities as Fig. [El but for the nu¬ 
merically simulated cluster halo g51. Original halo inclinations 
were set to i' = 40° and /' = 68°, and N\i = 5 iterations were 
used. The error curves are shown for the simulated halo with¬ 
out observational noise and using “radius-dependent smooth¬ 
ing”, and including observational noise and using the com¬ 
plete smoothing scheme. No noise was added to the emission- 
weighted temperature maps. The relative RMS Tew error was 
computed within a circle of radius 200 kpc around the map 
centre, while the density reconstruction errors were again deter¬ 
mined within the central 500kpc. Also for this numerical 
halo the errors are smallest for i « i'. 

For halos with an original inclination i' 90° or i' ^ 90° 
and for the analytic halos without noise, the minima of the er¬ 
ror curves for the emission-weighted temperature and the den¬ 
sity reconstructions are well dehned. Thus the quality of the 
reconstruction of such halos depends strongly on using the cor¬ 
rect inclination i = i' assumed in the deprojection. However, 
in such cases, the inclination angle is better constrained by 
the emission-weighted temperature maps. On the other hand, 
if the “symmetry” axis of the original halo is almost perpen¬ 
dicular to the line-of-sight, the minima of the error curves are 
usually broad, and finding the precise inclination i — i' for the 
reconstruction becomes less important. This can also be un- 
terstood from the fact that deviations are symmetric around 
i = 90°. For example, the halo with an inclination of /' = 68° 
shown in Fig. ini should exhibit minima at i = i' = 68° and at 
i = 180° — i'= 112° and a maximum in between. However, be¬ 
cause these three extremal points are close to each other, they 
start merging into one broad minimum. Note that the emission- 
weighted temperature maps can constrain the inclination angle 
in both cases to values where the errors of the reconstruction 
are close to their minima. 

The accuracy of inclination-angle estimates could be im¬ 
proved by using other independent information in addition to 
the temperature maps, such as data from weak and strong- 
lensing observations. Tensing observati ons allow reconstruc - 
tions of the lensing potential (see e.g. ICacciato et alJ boost , 
which is simply the suitably rescaled projection of the lens’ 
gravitational potential. A natural way to employ this for find¬ 
ing inclination angles is to assume hydrostatic equilibrium of 
the cluster gas in the potential well of the cluster and use the 
density and temperature reconstruction described above to ob¬ 
tain the gravitational potential of the reconstructed cluster halo. 


Its projection can then be compared to the lensing potential ob¬ 
tained from observations. Alternatively, one could use the de¬ 
projection algorithm discussed in Sect. l2.1l to obtain the three- 
dimensional gravitational potential from the lensing potential 
that was found from observations and compare it to the grav¬ 
itational potential corresponding to the density and tempera¬ 
ture reconstruction under the assumption of hydrostatic equi¬ 
librium. We shall explore these possibilities in a forthcoming 
study. 



0 10 20 30 40 50 60 70 80 90 

inclination i used for the reconstruction 


Fig. 12. Accuracy of emission-weighted temperature maps 
and densities of reconstructed analytic halos. The deprojec¬ 
tions started from maps obtained by projecting the analytic 
halo along a line-of-sight with inclination angles of i' = 40° 
and i' = 70°. Inclination angles i between 0° and 90° were used 
for the reconstruction. As expected, the best reconstructions are 
obtained for i « The errors are shown for deprojections from 
maps without observational noise and from smoothed maps 
with observational noise, and were averaged within a region 
of radius 500/;^* kpc around the map or halo centre. 


6. Summary and discussion 

We propose a new method for deprojecting hot gas 
in galaxy clusters which combines X-ray and thermal 
Sunyaev-Zel’dovich effect observations and reconstructs three- 
dimensional density and temperature distributions. We start 
from the iterative deprojection algorithm suggested by 
iRinnev et al ] ill99nli which employs Richardson-Lucy decon¬ 
volution and assumes axial symmetry of the physical quantity 
whose three-dimensional distribution shall be reconstructed 
from two-dimensional maps of its projection along the line- 
of-sight. 

This approach does not restrict the orientation of the sym¬ 
metry axis to be parallel to the line-of-sight, but the inclination 
angle between the symmetry axis and the line-of-sight is as¬ 
sumed to be known. This algorithm runs into problems when 
it is used to reconstruct strongly peaked distributions such as 
the X-ray emissivity of a galaxy cluster. There, one obtains 
spike-shaped artifacts through the centre of the reconstructed 
halo. We suppress the formation of such artifacts by introduc¬ 
ing a regularisation scheme for the iterative corrections used in 
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Fig. 13. Accuracy of emission weighted temperature Tew maps 
and densities for the reconstructed simulated halo g51. Maps 
obtained by projecting g51 along a line-of-sight with inclina¬ 
tion angles of i' = 40° and i' = 68° were used. The reconstruc¬ 
tion assumed inclination angles i between 0° and 90°. Best re¬ 
constructions are obtained for i « i'. The errors are shown for 
deprojections from maps without observational noise but using 
“radius-dependent smoothing”, and from maps with observa¬ 
tional noise on which the complete smoothing scheme was ap¬ 
plied. The RMS relative errors were obtained within a circle of 
radius 200 kpc around the map centre for the Tgw maps and 
inside a sphere of radius 500 kpc around the halo centre for 
the density reconstructions. 

the deprojection. Then, we generalise this algorithm to simul¬ 
taneously reconstruct three-dimensional distributions of sev¬ 
eral physical quantities by combining two-dimensional maps of 
projections which probe these three-dimensional distributions 
in different ways. 

Here, we discuss how three-dimensional density and tem¬ 
perature distributions of the intra-cluster medium can be 
reconstructed from combined X-ray and thermal Sunyaev- 
Zel’dovich effect observations. We test the method using syn¬ 
thetic data of analytically modeled and of numerically simu¬ 
lated galaxy clusters and discuss the quality of the reconstruc¬ 
tions and the impact of observational noise, cluster substruc¬ 
ture and deviations from axial symmetry. For numerical clus¬ 
ters which are of course not strictly axisymmetric, we use one 
of the principal inertial axes as the “symmetry” axis for the de¬ 
projection. 

Our main findings, if we neglect observational noise and 
assume that the inclination angle between the symmetry axis 
and the line-of-sight is known, are; 

- Spike-shaped artifacts of the deprojection are efficiently 
suppressed by the regularisation of the iterative corrections. 

- Densities and temperatures of the ICM of axisymmetric an¬ 
alytic clusters can be reconstructed very accurately from 
X-ray flux and Sunyaev-Zel’dovich effect maps. Errors are 
of the order of 1% unless the angle between the symmetry 
axis and the line-of-sight is small. 

- The three-dimensional density and temperature distribu¬ 
tions of hot gas in numerically simulated clusters, al¬ 
though not strictly axisymmetric, can still be reliably re¬ 


constructed. Relative errors reach roughly 10%. Smoothing 
of the X-ray flux and the Sunyaev-Zel’dovich effect maps 
can be used to suppress artifacts caused by subclumps. 

- Accurate gas density and temperature profiles can be ob¬ 
tained from the reconstructions. 

We then add photon noise corresponding to a total num¬ 
ber of 10''^ observed photons to the X-ray and observational 
noise corresponding to a four-hour ALMA Band 3 observa¬ 
tion to the Sunyaev-Zel’dovich effect maps, respectively. We 
smooth the maps before repeating the deprojections to suppress 
small-scale fluctuations which the algorithm would otherwise 
attempt to approximate and thereby reduce the quality of the 
reconstructions. From the repeated reconstructions of the an¬ 
alytic and numerical halos, this time including observational 
noise, we conclude: 

- Gas densities and temperatures of axisymmetric analytic 
halos can also be efficiently reconstructed from maps that 
contain observational noise. The relative errors of the re¬ 
constructions are about 5% to 10%. 

- The three-dimensional structure of the ICM of numerically 
simulated clusters can also be reliably reconstructed from 
X-ray flux and Sunyaev-Zel’dovich effect maps that contain 
observational noise. Relative errors reach roughly 15%. 

- Accurate profiles can be obtained from the reconstructions. 

- Five iterations are sufficient for the ICM deprojection. 
Using a larger number does not increase the quality of the 
reconstruction. 

For these deprojections, we assumed that the inclination an¬ 
gle between the symmetry axis and the line-of-sight is known 
beforehand. This will usually not be the case for observations 
of real clusters. In principle, one could try to And the inclina¬ 
tion angle by deprojecting the cluster assuming different values 
of the inclination angle and comparing the original X-ray and 
Sunyaev-Zel’dovich effect maps to those expected from the re¬ 
constructed halo. They should match best if the correct inclina¬ 
tion is used for the deprojection. Unfortunately, the minima of 
the deviations of the maps as a function of the assumed inclina¬ 
tion angle are quite broad and not always centred on the correct 
value. However additional data which are independent of the X- 
ray flux and the Sunyaev-Zel’dovich effect maps, can provide 
information on the inclination angle. We show that high-quality 
emission-weighted temperature maps which become more and 
more routinely available, can constrain the inclination angle of 
a cluster’s symmetry axis to values for which the quality of the 
reconstruction is close to its optimum. 
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